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\ Abstract 

Deutsch's algorithm for two qubits (one control qubit plus one auxiliary qubit) is extended 
Oh' to two d-dimensional quantum systems or qudits for the case in which d is equal to 2™, 

^ . n = 1, 2, . . . . This allows one to classify a certain oracle function by just one query, instead 

of the 2™ _1 + 1 queries required by classical means. The given algorithm for two qudits also 
solves efficiently the Bernstein- Vazirani problem. Entanglement does not occur at any step 
of the computation. 
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In the original Deutsch's problem the task is to ascertain a global property of an unknown 
but fixed Boolean function / : {0, 1} — ► {0, 1}, namely whether the function is constant (i.e. 
/(0)ffi/(l) = 0) or balanced (i.e. /(0)©/(l) = 1), where © denotes addition modulo 2. Suppose 
Cu we are given a "black box", or oracle, computing the function. Classically, it is obvious that 

one has to evaluate both /(0) and /(l) to solve the problem (that is, one needs to perform 
two function calls or queries). By allowing linear superpositions of single qubit states, however, 
Deutsch HQ devised a quantum algorithm which, with probability 50%, provides the correct 
answer for / (namely, whether it is constant or balanced) with a single evaluation of the function. 
. £^ \ An improved, deterministic version of Deutsch's algorithm has been given by Cleve et al. |2j 

which provides the exact answer for / in all cases. The speedup of the quantum algorithm 
over the classical one in this case is thus a factor of two. Deutsch's problem was generalized by 
Deutsch and Jozsa (DJ) P3 L to cover Boolean functions of the type / : {0, l} n — > {0, 1}. The task 
is again to determine whether / is constant or balanced, where balanced means that / yields the 
value for exactly half of the arguments and 1 for the rest. To solve this problem classically, it 
is necessary to get the function evaluated for 2 n_1 + 1 arguments in the worst case. Quantum 
mechanically, a single evaluation of / is sufficient [2j The exponential speedup achieved 
by the DJ algorithm stems primarily from the fact that it allows for quantum superpositions 
involving an exponentially large number 2 n of orthogonal states of n qubits. 

Another possible way to enhance the power of quantum computing is to increase the di- 
mensionality of the individual quantum systems involved in the computation. In this paper we 
provide an algorithm using two ci-dimensional quantum systems or qudits (one control qudit plus 
one auxiliary qudit) of dimension d = 2 n , which solves the generalized Deutsch's problem with 
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a single function evaluation. Note that the DJ algorithm employs n + 1 qubits (n control qubits 
plus one auxiliary qubit) to solve the same problem As both quantum algorithms (DJ's 

and ours) accomplish their goal with a single function evaluation we deduce that the computa- 
tional capacity of a single qudit of dimension d = 2 n is exactly the same as that corresponding 
to n qubits. Of course, this is directly related to the fact that the dimension of the state space 
of a single qudit with d = 2 n , is the same as the dimension of the state space associated with 
n qubits. The advantage of our algorithm resides in the fact that the control register consists 
of a single quantum system. Indeed it is important to emphasize from the start that, in our 
algorithm, entanglement does not take place at any stage of the computation. This contrasts 
with the DJ algorithm where, in general, entanglement does indeed occur between the qubits 
of the control register j^EUHj. Moreover, for a register consisting of a single ci-level quantum 
system, the steps of preparation, unitary evolution, and measurement involved in the quantum 
computation, should, at least for enough low values of d, be easier to perform than it would 
for the case of n distinct qubits. To substantiate this view, we must refer to previous work 
by Kessel and Ermakov 03 El HHj (see also Ref. JJ)- These authors, using the virtual-spin 
formalism 0IB1, showed how the four states of a single quantum particle (specifically, a spin-3/2 
nucleus) can be used to write and read two quantum bits of information, to prepare the initial 
state, and to implement a full set of two-qubit gates [3 El El- Analogously, they showed that 
a quantum particle with eight energy levels can store three qubits, and that three-qubit gates 
can be physically realized on one spin-7/2 particle [SJ^j. On the other hand, Ahn et al. ^2] 
experimentally investigated the storage and retrieval of information in the quantum phase of a 
coherent superposition state of energy levels in an A-state Rydberg atom, and, in particular, 
they demonstrated storage of numbers up to 2 N ~ 1 for N = 8. Furthermore, Ahn et al. 
suggested that a straightforward extrapolation of their results would allow numbers as large as 
2 100 to be stored in a single A-level atom, with A = 20. It should be noticed, however, that, 
although algorithms using single quantum particle as control register may not require entan- 
glement at any step of the computation, a physical realization of them would surely require an 
exponential overhead in some other recourse(s) for high values of d [H 1131 Hlj . 
Let us now define more precisely the problem under consideration. Let 

/:{0,l,...,d-l}-{0,l} (1) 

be a given function that maps each of the d arguments 0, 1, . . . , d — 1 to a one-bit value, with d 
being an even number 2 n , and n = 1, 2, . . . . The function is constrained to be either constant 
or balanced, i.e. it fulfills the property that 

/(o)©/(i)©...e/(d-i) = o, (2) 

for a constant function, or 

/(0)©/(l)©...e/(d-l)=d/2, (3) 

for a balanced one. Unless stated otherwise, in Eqs. (2) and (3), and in the remainder of this 
paper, © denotes addition modulo d. Given an oracle that evaluates the function for a given 
argument, our problem is again to decide, by queries to the oracle, whether / is constant or 
balanced. Note that the case of d = 2 corresponds to the problem originally considered by 
Deutsch pQ. Any (deterministic) classical algorithm for this problem would, in the worst case 
scenario, require d/2 + 1 function calls to know with certainty which kind of function we have at 
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Fig. 1: The network diagram for the generalized Deutsch's algorithm using two qudits. A 
final measurement result compatible with the control qudit being in the state |0) means 
that the function was constant. Any other output \x) (with x = 1,2, ... ,d — 1) means 
that the function was balanced. 



hand. As we shall presently see, there is a quantum algorithm involving two qudits that solves 
this problem with a single evaluation of /. This algorithm uses a quantum gate U / for two qudits 
that is a direct generalization of the /-controlled-NOT gate for two qubits used in the Deutsch 
algorithm. Let us denote the set of computational basis states of a qudit by {|0), |1), . . . , \d— 1)}. 
The operation of the two-qudit gate Uf is completely defined by its action on the computational 
basis for each qudit: 

\x)\y)^\x)\y®f{x)), (4) 

where \x) and \y) € {|0), \1) , . . . ,\d — 1)} denote the state of the control and auxiliary qudits, 
respectively. We may call the operation performed by the Uf gate a /-controlled-SHIFT opera- 
tion, for the effect of Uf on the auxiliary qudit when f(x) = 1 is to "shift" its state from \y) to 
the adjacent state \y © 1) (y = 0, 1, . . . , d — 1). 

The actual circuit implementing our algorithm is shown in Fig. 1. The initial state of the 
qudits in the quantum network is |0)|1), where the first (second) ket always refers to the state 
of the control (auxiliary) qudit. In the first step, both qudits undergo a unitary transformation 

(a qudit Hadamard gate) such that the state of the control qudit is transformed to 

|0)-|0H> = |0> + |l> + --- + |d-l>, (5) 

(apart from an unimportant normalization factor, which will be omitted in the following), 
whereas the state of the auxiliary qudit is transformed to 

|l)^|l H ) = |0)-|l) + -.. + |d-2)-|d-l), (6) 

with the + and — signs following each other alternately through the superposition. Then the 
state after the first two Hadamard gates is 

|0^>|1h>= J2(-l) y \x)\y). (7) 

x,y=0 

Next, we apply the unitary gate Uf on this state (the middle operation shown in Fig. 1). Of 
course, the actual computation of the function / comes from the action of this /-controlled- 
SHIFT gate. To determine the effect of Uf on the state (7), first note that, for each x = 
0, 1, . . . , d — 1, we have 

U f (X>imiy>) = Y J {-iy\x)\y®f{x)) = (-l)^\x)\l H ), (8) 
\y=o J y=o 
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since f(x) = or 1. The rightmost expression of (8) follows at once from the fact that, actually, 
the auxiliary qudit state \1h) is an eigenstate with eigenvalue (— l)^( x ) of the operator that sends 
the state \y) to \y © /(#))■ Therefore, the total state after the /-controlled-SHIFT is 



'd-\ 



£(-i) /(x) l*> |l ff > = |x/>|1h>. 



(9) 



\x=0 



From (9), we see that the state of the auxiliary qudit remains unchanged while each compo- 
nent \x) of the control qudit acquires a phase factor of (— l)^ x \ Clearly, for a constant function, 
the state \xf) is simply |0#). On the other hand, for a balanced function, the resulting state 
\Xf) will always be orthogonal to the state \0h), since now \xf) consists of an equally weighted 
superposition with exactly half of the |x)'s having a minus sign. Therefore, the two possibili- 
ties (namely, constant / or balanced /) can be reliably distinguished by means of a projective 
measurement P = \0h){0h\ on the control qudit. So, if / is constant (respectively, balanced), 
the probability of observing the control qudit in the state \0h) is 1 (0), so that a measurement 
filtering the state \0h) will with certainty give the result "1" ("0"). Alternatively, we may first 
apply a Hadamard transformation to the state \xf) in order to unitarily rotate the eigenbasis 
of the measurement into the computational basis {|0),|l),...,|d — 1)}. Thus, for constant /, 
the state of the control qudit after the last Hadamard gate becomes 1 0) , and for balanced / the 
resulting state becomes orthogonal to this, so that a subsequent measurement of the control 
qudit in the computational basis will distinguish these cases with certainty. (This is what is 
indicated in the last step of the diagram by the symbol M, Fig. 1.) As an example, for d = 4, 
the Hadamard operation has the possible matrix representation 
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Note that H\ = I4, and then, for example, we have that H^O) = |0) + |1) + |2) + |3) and, 
conversely, if 4 (|0) + |1> + |2> + |3» = |0). Similarly, we have that H A \V) = |0) - |1) + |2) - |3), 
and i?4 ( 1 0) — |1) + |2) — |3)) = |1). Analogously, for d = 8, the Hadamard transformation matrix 
has the possible form 
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Note again that H$ = 1$, so that Hg\x) = \xh) and, conversely, H%\xh) = for each 
x = 0, 1, . . . , 7. In general, the Hadamard operator acting on a single qudit of dimension d = 2™ 
is defined as 



d-l 



H d \x) = ±T(-ir x '\x'), 

^ d X >=0 



(12) 
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where 

X ■ X = ©j— q XiX^ (1^) 

denotes the bitwise inner product of x and x' expressed in the binary representation x = 
Y^2=o anc ^ x> = lLd=o x fi 1 ' w hh Xi,x[ £ {0, 1}. Please note that the symbol © in Eq. (13) 
denotes addition modulo 2. It is easily seen that, for each computational basis state \x), 
Hd(Hd\x)) = \x), so that the operator indeed fulfills Hj = Id- We further note that the 
matrix can be obtained from the matrix H2 = -j= ( } 1\ ) by performing n times the tensor 
product of i?2 by itself, i.e. i?2 n = ©"-#2- 

Provided with the general expression for the Hadamard transformation, Eq. (12), it is im- 
mediate to show that the quantum network in Fig. 1 can equally be used to solve the Bernstein- 
Vazirani (BV) problem ^1 (see also Ref. [2]): given an oracle which evaluates for some n-bit 
string a, the function 

f a (x) = a-x = ffi^ 1 a iXi (14) 

of domain x E {0, l} n , the task is to determine a by queries to the oracle. Note that there is 
a one-to-one correspondence between the n-bit domain {0, l} n and its decimal representation 
{0, 1, . . . , d — 1}, so that the function f a can equivalently be thought of as mapping the domain 
{0, 1, . . . , d— 1} to a one-bit value, as long as d = 2 n . Classically, to determine a it is necessary to 
perform at least n function evaluations, since a contains n bits of information and each classical 
function evaluation returns a single bit of information However, if the circuit in Fig. 1 is 
used with the function f a , then a can be determined with a single function evaluation. Indeed, 
in this case, and for the same initial state |0)|1) as above, Eq. (9) reads 

(j2(-l) aX \ x )) \1h) = \Xa)\lH). (15) 
\x=0 / 

Comparing the expression in parenthesis of (15) with (12), we can see that the control qudit 
state after the / a -controlled-SHIFT is simply the Hadamard transformation of the state \a), i.e. 
\Xa) = Hd\a). Therefore, after the action of the last Hadamard gate of the circuit in Fig. 1, the 
state \xa) becomes H^Xa) = H^\a) = \a). Thus, a final measurement of the control qudit in the 
computational basis will yield with certainty the integer a. 

Summing up, we have described a quantum algorithm using two qudits (one control qudit plus 
one auxiliary qudit) of dimension d = 2 n that solves with certainty the generalized Deutsch prob- 
lem, namely the problem of determining whether the function / : {0, 1, . . . , d— 1} — > {0, 1} is con- 
stant or balanced, with a single function evaluation followed by a measurement of the state of the 
control qudit. If we measure directly in the non-computational basis {|0#), . . . , \(d— 1)h)}, 
then finding the control qudit in the state \0h) unambiguously signals that the function was 
constant. Otherwise, if the control qudit is found in a state orthogonal to \0h), then / was bal- 
anced. On the other hand, if we measure in the computational basis {|0), |1), . . . , \d— 1)} (after 
applying a Hadamard transformation to the control qudit) and obtain the outcome correspond- 
ing to |0), then the function was constant. For any other output the function was balanced. 
Furthermore, we have shown that the same quantum algorithm using two qudits can be used to 
solve efficiently the BV problem. 

An important feature of our algorithm is that entanglement is never present throughout 
the quantum computation. Indeed, from Eq. (9), it is apparent that, after the action of the 
/-controlled-SHIFT gate, the control and auxiliary qudit states \xf) and \1h) are unentangled 
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from each other. Furthermore, each one of these states corresponds to a single quantum system 
and, therefore, entanglement is not an issue in the present algorithm. 1 Moreover, it should be 
noticed that our algorithm for two qudits can work equally well if the auxiliary qudit is replaced 
by a single qubit. In this case the gate Uf would correspond to a /-controlled-NOT gate defined 
by Eq. (4) (where now \x) £ {|0), |1), ... ,\d-l}} and \y) £ {|0), |1)}), and the qubit state \1 H ) 
would be simply |0) — |1), so that the joint (qudit plus qubit) state after the Uf gate would 
be given similarly by Eq. (9). Now, as the auxiliary qubit in state is not altered during 
the computation, we might further refine our algorithm by eliminating the auxiliary qubit while 
retaining the control qudit 5 . In this case, the function evaluation can be carried out by means 
of the /-controlled gate, whose action on the (computational) basis elements \x) of the control 
qudit is defined as 

U f \x) = (-l) f W\x). (16) 

Note that, when this refinement is made, we are indeed dealing with a single particle quantum 
computer! 

To conclude, it will be added that the logic behind the two-qudit gate Uf can be naturally 
extended to deal with multivalued functions of the type 

/:{0,l,...,d-l}-{0,l,...,d-l}, (17) 

in such a way that the operation of the U f gate is again defined by Eq. (4) , but now we let 
/ take the values f{x) = 0, 1, . . . , d — 1. It is readily seen that the effect of this generalized, 
multivalued gate on the state |0#)|1#) is given by 

c//|(M|1h) = |x/)|M, (is) 

where the transformed control qudit state \xf) is in turn given by 

\Xf)=^2(-l) f{x) \x). (19) 

x=0 

Again, Eqs. (18) and (19) are a consequence of the fact that \1h) is an eigenstate with eigenvalue 
(_l)/0) f th e operator that sends \y) to \y © f(x)), f(x) = 0, 1, . . . , d — 1. Now we define the 
following two classes of functions. We say that a function / has a constant parity when all values 
/(0), /(l), . . . , f(d — 1) are either even or odd. On the other hand, we say that / has a balanced 
parity when it yields an even value for exactly half of the arguments, and an odd value for 
the other half. A concrete example of a function with constant (balanced) parity for d = 8, is: 
(/(0), . . . , /(7)) = (4, 2, 0, 0, 0, 6, 2, 4) (respectively, (4, 2, 0, 0, 1, 1, 7, 5)). Then it is clear that if / 
has constant parity the state \\f) is simply |0#), whereas if / has balanced parity the state \\f) 
is orthogonal to |0#), since for this case exactly half of the terms in the superposition (19) have 
a minus sign. Thus, by suitably measuring the state of the control qudit, it can be determined 
with certainty whether / has constant or balanced parity, thereby showing that our algorithm 
using two qudits can also solve efficiently the generalized Deutsch's problem for multivalued 

x It is worth pointing out that the BV algorithm does not require entanglement even for the case where 
the control register consists of n distinct qubits [14| . Indeed, in this case, and for any n-bit string a, 
the expression in parenthesis of (15) can be written as a tensor product of n single qubit states. This 
factorization is a direct consequence of the above mentioned relation 7?2>» = <8>™i?2- This contrasts with 
the DJ algorithm for n control qubits, where entanglement is in general necessary 
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functions of the type (17). (Note that to solve this problem classically, d/2 + 1 function calls 
are required before determining the answer with certainty in the worst case scenario.) 2 

Finally, one might wonder whether the higher computational power of qudits could be ex- 
ploited to improve other existing quantum algorithms, or even to design entirely new ones. It is 
our hope that this will eventually be the case, and that the use of qudits and multivalued logic 
gates can serve as a valuable tool in the development of efficient quantum algorithms (see, in this 
respect, the work in Ref. ^H] where the authors derive a set of one- and two-qudit gates that are 
sufficient for universal multivalued computing, and show how such gates can be implemented by 
using d- level ions in the linear ion trap model). 
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2 For completeness, let us mention that by using a single 2™-level (respectively, 2 m ~level) quantum 
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function / : {0, 1, . . . , 2" - 1} -> {0, 1, . . . , 2 m - 1}, where n, m = 1, 2, . . . . 
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